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We study the possibility of creating a coherent population trapping (CPT) state, involving free 
atomic and ground molecular condensates, during the process of associating atomic condensate into 
molecular condensate. We generalize the Bogoliubov approach to this multi-component system and 
study the collective excitations of the CPT state in the homogeneous limit. We develop a set of 
analytical criteria based on the relationship among collisions involving atoms and ground molecules, 
which are found to strongly affect the stability properties of the CPT state, and use it to find the 
stability diagram and to systematically classify various instabilities in the long- wavelength limit. 

PACS numbers: 03.75.Mn, 05.30.Jp, 32.80.Qk 



INTRODUCTION 



The availability of atomic Bose-Einstein condensates 
(BECs) has made it possible to create, via photo- or 
magneto-association, molecular condensates. In photoas- 
sociation, a pair of free atoms is brought into a bound 
molecular state by absorption of a photon through a 
dipole transition. In magneto-association (or Feshbach 
resonance 0, ) , atom pairs in an open channel are con- 
verted into bound molecules in a closed channel through 
hyperfine spin interaction resonantly enhanced by mak- 
ing the energies of the two channels close to each other via 
the Zccman effect. Photoassociation creates molecules in 
excited electronic level, while magneto-association cre- 
ates molecules in high vibrational quantum state. In ei- 
ther case, the resulting molecules are energetically unsta- 
ble and suffer from large inelastic loss rate. Macroscopic 
coherence between atoms and molecules has been both 
studied theoretically 0, E H E ( also see a re- 
cent review article by Duine and Stoof |10|). and demon- 
strated experimentally |3(, nevertheless, long-lived stable 
molecular condensates have not been produced by direct 
association of atomic condensate • 

One possibility to overcome the difficulty is to employ 
the two-photon or Raman photoassociation model, as in 
the proposal forgenerating ground molecular conden- 
sates from atomic condensates through the 
stimulated Raman adiabatic passage (STIRAP) [IEG3- 
In spite of the nonlinear nature of the atom- molecule cou- 
pling in the pho toassociation process, it was shown by 
Mackie et al. [13| that in the collisionless limit, this model 
can support a coherent superposition of free atomic and 
ground molecular condensates. This superposition, be- 
ing immune to the inelastic processes associated with 
the excited molecular state, is called the nonlinear coher- 
ent population trapping (CPT) state, in direct analogy 
to the CPT state in the linear A three-level atomic sys- 
tem f]~8L ITg| . However, efficient operation of this scheme 




FIG. 1: The energy diagrams of three-level atom-molecule 
system involving free-bound-bound transitions. Conversion of 
atoms in \a) to quasi-bound molecules in \m) is accomplished 
by Feshbach resonance, while the coupling between \m) and 
the ground molecular state \g) is provided by laser light. 



typically demands very high laser intensities in order to 
overcome the weak coupling caused by the extremely 
small Franck-Condon overla p in tegral for the free atomic- 
bound molecular transition [l5j. 

To overcome the small Franck-Condon factor, a 

Feshbach-assisted STIRAP model is proposed [U l2l| . 
Here, conversion of atoms into quasi-bound molecules is 
accomplished by the much more efficient Feshbach pro- 
cess. To bring the quasi-bound molecules into a more 
deeply bound vibrational state, a two-photon Raman 
laser field is applied. In the limit where both optical fields 
are tuned far-off resonance from other electronically ex- 
cited levels, one can reduce the Raman coupling between 
the two molecular levels into an effective one-photon cou- 
pling. This leads to a three-level A atom-molecule system 
|22j| as shown in Fig.^ In this scheme, the possible source 
of decoherence is the decay of the quasi-bound molecu- 
lar state while the atomic state and the deeply bound 
molecular state are both stable. 

In the absence of particle collisions, we expect, from 
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the work of Mackie et al. [13J , that this system is capa- 
ble of supporting a CPT state, owing to the mathemati- 
cal equivalence to its photoassociation counterpart. The 
situation when the collisions are present is more com- 
plicated. It was generally believed that the nonlinear 
phase shifts arising from particle collisions would defeat 
the two-photon resonance condition |l3l Il5| . a prereq- 
uisite for the existence of CPT states and hence the key 
for the successful implementation of STIRAP. In a recent 
work |23|. we have shown that by appropriately chirp- 
ing the optical frequency and the Fcshbach detuning, the 
collision-induced nonlinear phase shifts can be dynami- 
cally compensated, hence generalizing the concept of the 
CPT state from the collisionless limit to situations 
where collisions cannot be ignored. In addition, we have 
also shown that collisions may lead to dynamical insta- 
bilities in the nonlinear CPT state, in contrast to the 
CPT states in the linear A-system which are always sta- 
ble. Avoiding dynamical unstable regimes is another key 
for successful implementation of the STIRAP process in 
nonlinear systems. 

The purpose of the current paper is to provide a sys- 
tematic study of the excitation spectrum of the gener- 
alized nonlinear CPT state, focusing on the character- 
ization of its stability phase diagram. To the best of 
our knowledge, a detailed stability analysis specific to 
the CPT state of the nonlinear A system has never been 
presented before. Our paper is organized as follows: 
In Sec. II, we present the Hamiltonian and the nonlin- 
ear CPT solution to the coupled Gross-Pitaevskii equa- 



tions. In Sec. Ill, collective excitations of the CPT state 
are calculated in the spirit of the Bogoliubov treatment 
[24l l25|. Sec. IV focuses on the stability analysis in the 
long-wavelength limit based on the excitation spectrum 
obtained in Sec. III. In Sec. V, we study the dynamics of 
a CPT state and show that dynamical instabilities may 
lead to self-pulsing of a condensate mode or growth of 
non-condensate modes. Finally, a summary is provided 
in Sec. VI. 



II. HAMILTONIAN, STATIONARY 
EQUATIONS, AND CPT STATE 

We choose to formulate our theoretical description for 
the Feshbach-based A-system as illustrated in Fig. ^ 
The exactly same description can be applied to the two- 
photon photoassociation-based A-system, due to the for- 
mal equivalence between the two models 0. In Fig.^ 
we use |a), |m), and \g) to denote the free atomic, quasi- 
bound molecular and ground molecular states, respec- 
tively. Levels |a) and |m) are coupled by a magnetic field 
through the Feshbach resonance with a coupling strength 
a and a detuning e (e is experimentally tunable via the 
magnetic field), while levels |ra) and \g) are coupled by a 
(effective) laser field with a Rabi frequency O and a de- 
tuning A. In second quantization, the total Hamiltonian, 
including collisions, takes the form 



H 



h J dr I E # « (~i v2 ) ^ (r) + \Xf Xij ^ (r) $ (r) & (r) ^ (r) 



+ e ft m (r) $ m (r) + | ( r ) ^ Q ( r ) ^ a ( r ) + ft. c .j + (A + e) $\ (r) ^ g (r) - - (r) $ g (r) + h.c 



(1) 



r 



where tpj (r) [ipi (r)] (i — a, m and g) is the creation 
(annihilation) operator of the bosonic field for species 
i at location r, and the terms proportional to \j rep- 
resent two-body collisions with = Anhai/rrii and 
Xij = Xji = 2irhaij / for i 7^ j characterizing the 
intra- and inter-state interaction strengths, respectively, 
(fli and ciij are s-wave scattering lengths, m t is the mass 
of species i with m m = m g = 2m a , and /Zy is the reduced 
mass between states i and j). We take, without loss of 
generality, both a and f2 to be real as their phase factors 
can be absorbed by a trivial global gauge transformation 
of the field operators. Here, we consider a uniform system 
and hence have dropped the external trapping potentials. 



Let us begin our formulation from the Bogoliubov ap- 
proximation, which amounts to decomposing ipi (r) as 



4>i (r) « ipi (r) + Sipi (r) . 



(2) 



where tpi (r) = (ipi (r) > is the condensate wave function 

and Stpi (r) is a small fluctuation field operator that obeys 
the usual bosonic commutation relation 



#i(r),M(r') =5i^(r-r') 



(3) 



As a first step in the Bogoliubov approach, we put Eq. 
into the grand canonical Hamiltonian lfT)l 



H = H - fyi dr 



n 



2^4 



(4) 
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where the chemical potential hfi is introduced to conserve 
the total particle number, and expand Eq. Q up to the 
second order in fluctuation operators as 



(5) 



where the subscript represents the order in Sijii (r) . 

The zeroth order term, H^ ', is a constant depending 
on the wave functions ipi and is irrelevant to our interest. 
The first-order term, H.^ 1 ', is required to vanish in the 
Bogoliubov formalism. This requirement leads to a set 
of coupled stationary equations 



A Qi \tpi\ J ip a + aijjm^*, (6a) 




(6b) 



]>> g4 |^| 2 V ff + (A + e)V 9 



(6c) 



which represent the generalized time-independent Gross- 
Pitaevskii equations (GPEs). Finally, the second-order 
term, W 2 ', which has a more complicated structure, will 
be studied separately in the next section where we cal- 
culate the spectra of the collective excitations. 

We now turn our attention to the homogeneous (zero- 
momentum) solution to the GPEs. To proceed, we sep- 
arate each wave function into an amplitude and a phase 
according to 



(7) 



To seek the CPT solution, we take ip m — 0. After sub- 
stituting Eq. J7J) into Eq. ©, we find that under the 
conditions 



a \ip a 



29 a , 



(8) 



Eq. © possesses a solution consistent with the require- 
ment that ij) m — 0. The explicit form of this solution can 
be readily found by combining Eq. (JSJ ) with the particle 
number conservation 



where n is the total particle density, to yield the particle 
density distribution as 



m 2 = o, 
k\ 2 = - 



2)i 



l+\ 1+8 



V " ) 



ot\/n 



(10a) 
(10b) 



(10c) 



1 + 1/1 + 8 



A check for consistency leads to the chemical potential 

M = A Q |^| 2 + A ag |^| 2 , (11) 
where Ai = Am, and the laser detuning 



(2A afl 



A g )U°| 2 + (2A a -A ag )|^| 2 , (12) 



\^ a \ 2 + 2\%\ 2 =n, 



(9) 



which is necessary in order for Eqs. i|10[l to be a stationary 
solution of Eqs. ©. 

The vanishing of the population in the quasi-bound 
molecular state indicates a destructive interference in the 
excitation to this state. A system prepared in the state 
represented by Eqs. (fTU|) is not subject to the particle 
loss suffered by the quasi-bound molecular state. This 
situation is reminiscent of the CPT state in a linear A- 
type atomic system where the atoms are immune to the 
spontaneous emission 0,0]. For this reason, we regard 
the steady state represented by Eqs. (|1(J|) as the non- 
linear matter-wave analog of the CPT state in a linear 
A-system. 

Interestingly, just as in the linear case, the CPT so- 
lution described by Eqs. I|10|) depends only explicitly on 
the two coupling strengths, a and SI, not on the colli- 
sional parameters. The latter, however, play important 
roles as evidenced by Eq. (|12fl . Note that in the absence 
of collisions, Eq. i|12fl degenerates into A + e = 0, which 
is just the ordinary "two-photon" resonance condition. 
Thus, Eq. I|12|l can be regarded as the generalization of 
the two-photon resonance in which the nonlinear phase 
shifts due to particle collisions have been compensated. 
This CPT solution is therefore a nonlinear generalization 
of the one found in the collisionlcss limit 131. 



III. COLLECTIVE EXCITATIONS OF THE CPT 
STATE 

In this section, we want to calculate the excitation 
spectra of the CPT state, especially, the spectra in the 
long- wavelength limit, since they determine many impor- 
tant properties, including the stability, of the CPT state 
at low temperature. 
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A. Bogoliubov Equations 



where we have defined 



To determine the equations for the collective excita- 
tions, we first take advantage of our system being homo- 
geneous and move to the momentum space through the 
expansion 



Sipi (r) 



J— fdpj^Ciip), (13) 
7m) J 



(2nh) 



where (p) [cj (p)] is the annihilation (creation) opera- 
tor for a particle of species i with momentum p and obey 
the standard bosonic commutation relation 



Ci(p),ct(p') =Sij6 (p-p') 



(14) 



in accordance with Eq. J2J). The third term Ti.^ in 
Eq. (JSJ for our CPT state can now be put into a compact 
form 

nm 

+ [c n (p')B nm c m (-p')+h.c.}}, (15) 

where A and B are both real and symmetric matrices, 
given by 



A = 



B 



a\ip° a \ e f + e' -§ 



A Q |<| 2 A Q9 |^||<| 


A afl |v2|KI x 9 \<\ 2 



(16a) 



(16b) 



e p = p 2 /(2m a h) , 

e' = £ -A a |^| 2 -A g |^| 2 , 



A„ 



2A a ^ ^mg ■ 



(17a) 

(17b) 
(17c) 
(17d) 



The process of finding the quasiparticle spectrum be- 
gins by changing Eq. I|15fl into a diagonalized form 

7Y (2) = ft J dp^uji (p) b\ (p) k (p) + c-number, (18) 



with the generalized Bogoliubov transformation 

k (p) = u » (?) £ J (P) + H u y O 5 ) £ ] (~P) ' ( 19 ) 



where LOi (p) is the quasiparticle frequency, Uij (p) and 
(p) are the transformation coefficients, and bi (p) and 
frj (p) are the respective annihilation and creation oper- 
ators for a quasiparticle of frequency uji (p) and momen- 
tum p. To obtain the equations for w$ (p), ity (p) and 
Wy- (p), we first note that 



h{p),n (2) ] = ^ (p) k ( P ) . 



(20) 



Inserting Eq. (|T§|| into Eq. l(2U|) and with the help of 
Eqs. i|14|) and (|15fl along with the symmetric properties 
of A and B, we arrive at 



^2 X! { 1^™™"™ ~ B mn v m ] c m (p) + [B rnn u ln - A mn v ln ] & m (-p)} = LU l ^ [u im c m (p) + u im c^ (-p)] . (21) 



Equating the coefficients of c m (p) and & m (— p) on both 
sides of the above equation, we arrive at a 6x6 matrix 
equation 



where the coefficients are given by 



A -B 
B -A 



Ui (p) 



(22) 



where 



l^i (^ia ) l^im ? ^ig ) ? — (^ia; ^im^ ^ig] 

A simple manipulation can show that uif(j>) are the 
eigenvalues of the 3x3 matrix (A + B)(A — B). hence 
(p) satisfy a cubic equation with the form 



(w?) - oi (w, 2 ) + a 2 



«3 = 0, 



(23) 



Ol 

O.'i 



= oio + ane p + ai 2 e p , (24a) 
= »20 + a 2 ie p + a 22 £p + a 23 £p + a 2 4£p , (24b) 



a 3 i£p 



i 



a 32 e p 

,6 



a 33 e„ + a 34 e 



(24c) 



The coefficients depend on the coupling and the col- 
lisional constants and are in general quite complicated. 
We only list three of them below 
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/2 1 
= £ + 2 


a 2 o 


= -Iftt 2 
16 V 


«:S1 


= - m 2 



-(fi 2 + 4a 2 |^ a | 
^ (f> 2 + 4a 2 |^| 2 ) 2 - 0.5e'a 2 |V>°| 4 (A s - 4A ag + 4A a ) , 



) x [A a fi 2 +4a 2 (\ g \^ g \ 2 + KM 2 ) -8e'(A Q A 9 -A 2 g )|^| 



(25a) 
(25b) 
(25c) 



as they will be explicitly referenced in later discussions. mentum with 



B. Excitation Spectra 

Equation (|23(l has a cubic form with respect to uJf and 
hence can be solved analytically. However the analyti- 
cal expressions are in general too complicated to provide 
much physical insights. Since we are most interested in 
the long-wavelength excitations, we proceed to take a 
perturbative approach by expanding u>i(j>) in the limit 
of small p, or equivalently, small ^fe^. The perturbative 
solution thus takes the form: 



{p) 



,(o) 



(26) 



In order to bring out the effects of collisions more clearly, 
let us first consider the collicionlcss limit. 



1. Excitation Spectra in the Collisionless Limit 

In the absence of collisions, the matrix B becomes null. 
The excitation frequencies Wj(p) are eigenvalues of the 
Hermitian matrix A and therefore must be real. This 
means that there is no dynamical instability in spite of 
the nonlinear nature of the coupling between the atomic 
and the quasi-bound molecular state. 



To the zeroth order, we find that ui 
tions H| 



(0) 



has three solu- 



(0) 



aw ± 



v /a 2 -4a 20 ) /2, (27) 



where the subscripts (0,—,+) are used to denote the 
three branches of the collective modes. Clearly the 0- 
branch is gapless, while the (±)-branches are gapped. 
The reality of uj^ is guaranteed because here both aio 
and a2o are positive with a 2 > 4a2o [see Eqs. (|25a|) and 



(2) 



7(2) 



» 2 +2« 2 i^r 

2 +4a 2 |<| 2 



± 







2 ( W f ) 2 ] 


2 


(4 0) ) 2 + 


n 2 

4 


+ « 2 !€| 2 ' 



(28a) 



(28b) 




FIG. 2: The quasiparticle dispersion spectra for n = 1, fi = 1, 
and e = —2. Solid lines represent the collisionless case while 
dotted lines are for a system with A a = 0.625, A m = A = 



Aam — Aag — 



0.3A a . Units are defined in Sec. II VI 



(125 bfl . in the absence of collisions.] 

For all three branches, we have found that d. 
Hence the next leading order in Ui (p) is quadratic in mo- ification of the gapless 0-branch. While d±' still re- 



(i) 



0. 



Solid lines in Fig. [^represent the three branches in the 
collisionless limit. They are obtained by directly solving 
the cubic equation l|23|l . and are in good agreement with 
the perturbative results (not shown). 



2. Excitation Spectra with Collisions 

The spectra including the collisional terms can be sim- 
ilarly obtained. The zeroth order solutions, o^ "*, have 
exactly the same forms as in the collisionless limit given 
by Eq. I|27|l . Therefore the 0-branch continues to be gap- 
less. Due to the collision-modified coefficients , the 

two zeroth order gapped modes are now, however, no 
longer guaranteed to be real. Complex mode frequencies 
lead to the dynamical instability, which we will discuss 
in detail in the next section. 

Another important effect of the collisions is the mod- 

(i) 
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mains to be zero (hence the two gapped branches are 
still quadratic in p to the leading order), now takes 
a finite value given by 



7(1) 




(29) 



Hence the gapless branch, in the long- wavelength limit, 
is linear in p and represents the phonon mode with the 
speed of sound determined by Eq. (|29|l . The dotted lines 
in Fig.[21represent the spectra with collisions in the stable 
regime. Once again, they are calculated by solving the 
cubic Eq. I|23[l , but are found in good agreement with the 
perturbative results (not shown). 

According to Eq. I|29l) , whenever 031 and 020 are oppo- 
site in sign, luq becomes imaginary for finite p in the long 
wavelength limit, indicating an unstable 0-branch. We 
now turn to a detailed discussion of the collision-induced 
dynamical instability. Note that from Eq. I|29|) . it may 
seem that ojq is singular for any finite momenta when 

a20 = 0. However, when a2o = 0, lu^ and uj'f ' are de- 
generate according to Eq. (|2*T|) . Thus, this singularity is 
only an artefact, indicating the failure of the nondegen- 
erate perturbative calculation. Indeed, exact solutions 
of Eq. I|23[) do not exhibit such a singularity. Despite of 
this complication, the perturbative results can be used 
to identify the unstable regimes quite accurately. 



IV. 



CLASSIFICATION OF DYNAMICAL 
INSTABILITY 



Studies in Sec. IIII B~2l suggest that complex excitation 
frequencies and hence dynamical instabilities can occur 
under the following situations: 

1. When 031 and 0,20 possess opposite signs, 0-branch 
becomes unstable [see Eq. (|2^|) ]; 

2. When a 2 o < 0, w_ becomes complex [see Eq. <|2*7|) ]: 

3. When a 2 < 4ft2o, both (±)-branches become un- 
stable [see Eq. l|27|)] . 

Before moving toward, let us clarify our unit system. 
We consider systems with particle densities in the or- 
der of n = 5 x 10 20 m~ 3 , and a fixed to 4.22 x 1CT 6 
m 3/2 g -i^ which corresponds to the atom- molecule cou- 
pling strength for the 23 Na Feshbach resonance at a mag- 
netic field strength of 85.3 mT 27). We then adopt a 
unit system in which no is the unit for density, a/ \Zn~o = 
1.88724 x 10 _16 m 3 s _1 the unit for collisional parameters 
(Ai, Ay), a^/no — 9.436 x 10 4 s _1 the unit for frequen- 



cies (w,fi,e), and y (2mhay/no) = 8.72 x 10 28 kg m/s 
the unit for momentum p. Furthermore, for simplicity, 
we only consider the situation where all the scattering 
lengths are positive. 



Our goal is to identify the unstable regimes in the pa- 
rameter space spanned by e and O, while the laser de- 
tuning A is always fixed by the resonance condition rep- 
resented by Eq. i|12|) . For cases 1 and 2, we can identify 
two threshold values for the Feshbach detuning e h and 
e\ , given by 



th 





2 


KA 3 






2 


fa 9 





( Q 2 



1 



4A Q +Xn- 4A og 



x a n 2 



+ 4(A ff |<| 2 + A ag |V> a °| 2 ) 



8|^°| (A a A fl 



Kg) 



(30) 



(31) 



which are obtained from the conditions 020 = and 031 
0. respectively. Case 3 will be discussed at the end of this 
section. 



A„=A 




FIG. 3: Division of the positive quarter of (A s , A a )-space into 
three regions by a straight line: A fl = — 4A a + 4A ag and a 



hyperbola: A a A s 



A 



It is not difficult to see that the signs of the denomina- 
tors of the last terms of l{3"0"|l and l|3"T|) . 4A a + A g — 4A ag and 
A a A 9 — A 2 ff , respectively, will play crucial roles in locating 
the instability parameter regime. For this reason, we first 
decompose the (A s ,A a )-space into three regions divided 



by a straight line, A„ 



-4A a + 4A a „, and a hyperbola, 



A a A 5 = A 2 ff , as shown in Fig. The line is tangent to 
the hyperbola at the point (A g = 2\ ag ,X a = 0.5A Qg ). The 
three regions are defined as follows: Both 4A a + A g — 4A ag 
and A Q A g — A 2 g are negative in Region I and both are pos- 
itive in Region II, while in Region III, 4A a + A g — 4A ag < 
and A a A s - \ 2 ag > 0. 

Next, we choose to construct, for each region defined 
above, two instability maps in the (f2, e)-space; one for 
the instability due to the imaginary sound speed (case 
1), and the other for the instability due to the com- 
plex gap (case 2). We represent them by the shaded 
areas of the diagrams in the first and second columns 
of Fig. 21 respectively. Each set of diagrams in row 1, 
2 and 3 are built upon the collisional parameters corre- 
sponding to Regions I, II, and III of (A g ,A a )-space, re- 
spectively. We construct these maps based on two prin- 
ciples. First, it can be shown from Eqs. I|30|) and l|31|) 
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Case 1 



Case 2 



III 




FIG. 4: Classification of the instabilities in the limit of low 
momentum. The shaded areas of the diagrams in column one 
and two correspond to the instabilities due to the imaginary 
sound speed and the complex gaps, respectively. Each set 
of diagrams in row 1, 2 and 3 are the instability maps for 
systems operating in region I, II, and III as specified in Fig. 
3, respectively. 



that ef - e ' 1 oc 1/ (4A a + X g - 4A ag ) (A a A 9 - \ 2 ag ) with 
a positive definite proportionality. As a result, we have 
ef > e tjh in Regions I and II, and e\ h < Eq in Region 
III. Second, according to Eq. (jSTJ [Eq. (jSU)], 4 h [ef] is 
bent downward if 4A a + X g — 4X ag < [A a A g — A 2 g < 0] 
and vice versa. Note that the instability of case 2 is solely 
determined by eg' 1 ; the presence of e[ h in the second col- 
umn of Fig. 0] is to make the comparison of the two maps 
easier. For example, one can easily identify that the two 
instability maps in region II overlap in the band between 



eg 1 and e\ h . 

Let us now consider a specific atom-molecule sys- 
tem composed of 23 Na. The atomic s-wave scattering 
length for sodium is around a a = 3.4 nm psj which 
yields A a = 0.625. No precise knowledge about the 
scattering lengths involving molecules exists. Assum- 
ing they are on the same order of magnitude as a a , 



we take A„, = 



K = 



X„ 



= A 



my 



— X a „ — 0.3A Q . 



This set of parameters puts the system in Region II of 



Ah 



and 



Hi 



Fig. |U According to Eqs. ljA3j) and (|A"4l 

approach 0.5A g + 2a 2 / (4A a + X g - 4A ag ) = 1.126 and 
0.5A g + 0.5a 2 A g / (A a A g - A 2 g ) = 1.237, respectively, in 
the limit of small tl. 

A search using Eq. (|23(1 above these thresholds and 
at small momenta indeed leads to two branches whose 
frequencies have finite imaginary parts for certain val- 
ues of fi, as illustrated in Fig. OJa) and (b), respectively. 
Each curve in Fig.E^a) is restricted between two thresh- 
olds, both of which increase as e, evidently resembling the 



0.02 



0.015 



S o.oi 



0.005 







(a) 


J 









0.5 



1 

Q 



1.5 




FIG. 5: Im[u; (p)] as a function of Q for a fixed momentum 
p = 0.005. The curves from left to right correspond to e = 1.3, 
e — 1.5, e = 2, and e = 3, respectively. Other parameters are 
n = 1, X a — 0.625, and A m = 
0.3A a . Units are defined in Sec. II VI 



A 9 = 



A„ 
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Xa 



instability originating from the imaginary sound speed 
(case 1). Thus, in the limit of small p, we expect it to ap- 
proach d\pp. In fact, one can trace the sharp asymmetry 
of each curve to the fact that according to Eq. (|2*§|l . cfip 
vanishes at the left threshold where 031 is close to zero, 
while becomes singular at the right threshold where 020 
is close to zero. In contrast, each curve in Fig. EJb) has 
only one threshold, which matches quite closely to the 
right threshold of the corresponding curve in Fig. p [a) . 
This along with the fact that the value of the mode fre- 
quency is rather insensitive with respect to p (as long as p 
is not too large) indicates that it represents the instabil- 
ity from the complex gapped branch of case 2. Therefore 
these exact calculations are in good agreement with the 
insights gained from the instability map shown in Fig.0] 

In order to have a more detailed look at various insta- 
bilities, we make a three-dimensional plot in Fig. for 
the imaginary parts of the two (non- vanishing) branches 
of excitation frequencies. A glance at this 3D view fur- 
ther confirms the instability map of Fig. 4: one can easily 
associate the large bump in Fig. Ela) to the instability 
of case 2, and the curved band in Fig. HJb) to the in- 
stability of case 1. The shapes of these unstable regions 
completely match with those shown in Fig. 0] 

However, Fig.^Jb) also reveals a new feature consisting 
of a narrow strip along the £1 dimension. We can trace 
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FIG. 6: (Color online) A three-dimensional view of Im[u} (p)] 
as a function of SI and e for two branches of roots to Eg, 1231 . 
Parameters are same as in Fig. |K| Units are defined in Sec. 



this new feature to the source of instability associated 
with case 3 which has so far been ignored. An analysis 
with the help of Eqs. (|25a[l and (|25b|) indicates that the 
condition 020 > Qio/4 amounts to 

e' 3 + (Sl 2 + 4a 2 |^| 2 ) e' 

+2a 2 |^| 4 (A s - 4A ag + 4A a )] < 0, (32) 

from which we find that to a good approximation, the 
unstable e is bounded between 



ei 

£2 



A a |^| 2 + AJ^| 



2a 2 |^| 4 A 9 - 4A ^ + 4A - 



n 2 + 4a 2 \^\ 



(33a) 
(33b) 



It can be shown that e± and £2 are not very sensitive 
to f2, and the they merge together in the limit of either 
small or large f2, which explains the shape of this new 
instability region. 



V. DYNAMICAL SIGNATURES OF STABLE 
AND UNSTABLE CPT STATES 

In order to determine the fate of an unstable state, 
we must examine its dynamical response to fluctuations, 
which can be simulated using the time-dependent version 
of the GPEs, which are obtained by replacing the left 
hand sides of Eqs. JfJJl by iipi. 



We take a CPT state associated with a Rabi fre- 
quency flo and Feshbach detuning e to be the initial state. 
To simulate small fluctuations, we modify the Rabi fre- 
quency as 



fi = fl + [5tt + 2(50i sin (px/H)] , 



(34) 



where SQq and <5f2i are small perturbation amplitudes. 
Note that a position-dependent perturbation represented 
by the Sfli term can cause spatial deformation of the uni- 
form CPT state, making it possible to trigger the insta- 
bilities due to perturbations of finite momentum. 
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FIG. 7: The time evolution of the particle number densities in 
the condensate mode probed by a zero-momentum perturba- 
tion with 5Q = 0.001 for (a) fi = 1.1 and (b) Q = 0.9, while 
the Feshbach detuning is fixed at e — 3. Other parameters 
are the same as in Fig.|S] Units are defined in Sec. II VI 

To test the instability of case 2, we take the same colli- 
sional parameters as in Fig-EJ but fix the Feshbach detun- 
ing to e = 3. We apply a zero- momentum perturbation 
with an amplitude <5S1 = 0.001 and <5Qi = to the CPT 
state with (a) flo = 1.1 and (b) Oo = 0.9, respectively. 
By inspecting the e — 3 curve of Fig.[5fb), which survives 
even when p = 0, we anticipate that the dynamics of the 
system will change from a stable to an unstable one as f^o 
varies from 1.1 to 0.9. The results, shown in Fig.[7{a) and 
(b), clearly support our analysis. The stable CPT state 
does not respond significantly to the small perturbation; 
while Fig. EJb) shows a large-amplitude self-pulsing be- 
havior, characteristic of an unstable system. In addition, 
there is a noticeable population in the unstable molecular 
state. Such self-pulsing oscillations are also familiar fea- 
tures in the study of unstable nonlinear optical systems 
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To demonstrate the instability of case 1, we need to 
apply a finite- momentum perturbation since for case 1, 
the imaginary part of the excitation frequency loq van- 
ishes at p = 0. The finite- momentum perturbation will 
couple together different momentum modes. To simulate 
this process, we adopt the Floquet technique and expand 
the wave functions as 

-f-oo 

A(r,t)= J2 ^{l,t)^ x/h , (35) 

/— — oo 

where \<f>i (l,t)\ 2 represent the particle density of the ith 
I 

. dfa (n) 

I : 

dt 

. dfan (n) 

l : 

dt 



. d<p g (n) 
1 dt 



The initial condition is then equivalent to 

M0,0)= |VS|.^m (0,0) =0,^(0,0) = |^| , (38) 

and <f)i (I ^0,0) — 0, where |?/>°| and \?pg\ are deter- 
mined from Eqs. i|l(J I with f2 being replaced by fio- 
We now switch to a set of collisional parameters with 

A a = 0.625, X m — A ff — A m g = 0.5A a , \am = A a , 

and X ag = 1.4A a , which puts the system in Region I, 
where the two instabilities regions of case 1 and 2 do not 
overlap in the same parameter space (see the two dia- 
grams in the first row of Fig. 0J. We take e = —4 and 
J2q = 1, which puts the system in the unstable region 
of case 1. This allows us to demonstrate the instability 
due solely to the finite-momentum perturbation. Indeed, 
simulations show that the system is stable against zero- 
momentum perturbations but is unstable when finite- 
momentum perturbations are applied. A dynamical re- 
sult of a finite-momentum perturbation is displayed in 
Fig. Q3] which shows that the instability leads to an irre- 
versible population transfer from the condensate mode 
(p = 0) to non-condensate modes (p ^ 0). The initial 
growth of the non-condensate modes is approximately 
exponential but quickly becomes rather chaotic, which is 
typical for an unstable multi-mode system |30| . 



species with momentum Ip. Here, I is an integer and 
p = 2irh/L, with L being the length of the condensate 
along x dimension. Then the orthonormality condition 
for the momentum modes reads 



J e i{ l - l '>*/ h dx = 6i tV . (36) 



This enables us to transform the dynamical GPEs for 
ipi (r,t) into the coupled modal equations for fa (l,t) 



(37a) 
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FIG. 8: The time evolution of \fa (I, t)\ 2 (I = 0, 1, and 2) after 
a finite-momentum perturbation with Sfli = 0.001 at p = 0.01 
[5Q.0 = 0) is applied to a system initially prepared in a CPT 
state of e = —4 and flo = 1- We include modes ranging 
from I = — 10 to 10 in the simulation. Other parameters are 

A a = 0.625, A m — Xg = \mg — 0.5A Q , A Q m = A a , \ a g = 1.4A a , 

and n — 1. Units are defined in Sec. II VI 

VI. CONCLUSION 

To summarize, using a generalized Bogoliubov ap- 
proach, we have calculated the collective excitation spec- 
tra of the CPT states involving atomic and stable molec- 
ular condensates. We have found, through our analysis 
of the excitation spectra, that collisions involving atoms 



-e p n 2 fa (n) + ^2Yl Xas< ^ s (0 # (0 0a{n + l' -I) + a^(j) m (I) fa a (I - n) , 



LI' s 



-yn 2 m (n) + X rns<Ps (0 4>* s (0 <t>rn (n + I' - I) + efan (n) 

1,1' s 



+9 J2^ a ( n ^) 



(n)-^[fa (n-l)-faj (n + 1)] , 



■ e fn 2 fa (n) + Y,Y, X ^ M # (0 ^ (n + l' — l) + (A + e) fa (n) 

1,1' s 



n + 6n sn x 



(n-l)-ct> m (n + l)] 
I 
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and stable molecules strongly affect the stability proper- 
ties of the CPT state. We have developed a set of analyt- 
ical criteria for classification of various instabilities in the 
long- wavelength limit. We have shown that these criteria 
can be used to accurately identify the unstable parameter 
regimes. In this paper, we studied a homogeneous sys- 
tem. It will be desirable in the future to extend this work 
to the inhomogeneous case where trapping potentials are 
present, as in typical experimental situations. Neverthe- 
less, our work here should provide important analytical 
insights into the much more complicated inhomogeneous 
problem, for which the only viable treatment is perhaps 
through numerical simulations. 

Just as the usual CPT states in linear three-level A- 
systems are behind many important applications in quan- 
tum and atom optics [3l|, we anticipate that the nonlin- 
ear CPT states will also play important roles in appli- 
cations where the coherence between atom and molecule 
condensates becomes critical. Evidently, the success of 
these applications depends on our ability to avoid the 
unstable regimes of the CPT state. For example, high ef- 
ficiency coherent population transfer between the atomic 
and molecular condensate using STIRAP can only be 
achieved if these unstable regimes are not encountered 
during the transfer process |23j . With this in mind, we 
believe that our work will have important implications 
in the coupled atom-molecule condensates |7j, a system 
of significant current interest, with ramifications in other 



nonlinear systems. 
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APPENDIX A: ASYMPTOTES 

A qualitative understanding of the instability map 
shown in Fig. 4 depends on the asymptotic behavior of 
the threshold Feshbach detunings and e\ h which we 
show below. In the limit of large O/a 

^ " Aa+ 8(4A a + A s -4A QS ) (S) ' (A1) 

«5» » A a+ sfA f (A2) 
8(A a A 9 -A 2 g) \aj 

and in the limit of small Q/a 



Ah 



Ml 



0.5A o 



0.5A o 



2a 2 



4A a + X g - 
0.5a 2 X g 

X a X g — X a 



4A„ 



2A a - A fl 
2V2 



V2a 2 



4A a + X Q 



4A„ 



(2Aqg - Ag) + AgAg (Ag 

2V2(A 2 g-A a A s 



2A Q ) - 2a 2 A ag fl 



(A3) 
(A4) 
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